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ABSTRACT 

A study of correlations in tractable multiparticle cascade models in terms of 
wavelets reveals many promising features. The selfsimilar construction of the 
wavelet basis functions and their multiscale localization properties provide a new 
approach to the statistical analysis and analytical control of hierarchically or- 
ganized branching processes. The exact analytical solution of several discrete 
models shows that the wavelet transformation supresses redundancy in the cor- 
relation information. Wavelet correlations can be naturally interpreted as corre- 
lations between structures (clumps) living on different scales. 
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1 Introduction 



Multiparticle production is often dubbed in terms of branching processes. In this context, 
selfsimilarity and scaling are particularly appealing concepts since they constitute an in- 
triguingly simple and universal organisation principle of the underlying dynamics. They are 
found in a large variety of complex phenomena. 

Simple discrete hierarchical branching models with selfsimilar density fluctuations were 
first introduced in the description of energy dissipation in fully developed turbulence 
They combine an independent hierarchical evolution of different branches with a random 
multiplicative nature of density fluctuations. In multiparticle processes such selfsimilar cas- 
cade models serve as paradigmatic toy models for intermittency Intermittency is now 
experimentally well established in soft multiparticle interactions and is now being explored 
in hard (perturbative) multiparton interactions. 

Although these branching models describe local density fluctuations, they are convention- 
ally analyzed (globally) in integrated form as is the case for the experimental data: The 
underlying concepts of scaling and selfsimilar fluctuations are revealed in a multifractal (mo- 
ment) analysis f|, which globally averages over the differential (local) correlation structure. 
A fully differential correlation analysis itself would be a more fundamental approach to re- 
veal selfsimilarity and scaling and deviations from these. However, correlation functions of 
higher than second order are difficult to measure and to visualize. Nevertheless, we prefer to 
stick to a fully local description of hierarchical branching models and their analysis in terms 
of differential correlation densities. 

One word about the branching models used in the course of this paper: Alltogether they 
serve as (oversimplified) discrete approximations to realistic QCD parton cascades in the 
perturbative regime. We do not intend to overstress this relation, but rather prefer to 
study new correlation techniques demonstrated by relatively simple branching models. 
We restrict the discussion to one-dimensional discrete branching processes only. Higher 
dimensional branching processes can be treated in the same way as the one-dimensional 
case and do not change the conclusions drawn. Discreteness which is introduced into the 
branching process by allowing branching only for discrete "time" steps of the cascade, is also 
not a fundamental shortcoming compared to continuous processes as long as the step width 
stays small. We are aware that a complete continuous formulation of the branching process 
needs a functional approach. 

In this paper we essentially convey three messages. First, we generalize the univariate 
(moment) generating function to a multivariate (correlation density) generating function 
and derive its evolution equation. Message two unveils a secret unnoticed so far in moment 
(multifractal) analysis, but which is essential for comparison of model predictions with data: 
For cascade models with global density fluctuations on top of local density fluctuations, a 
careful distinction has to be made between moments obtained from a (theoretical) forward 
evolution of the cascade and moments obtained from an (experimental) backward analysis. 
Last not least, as message three we propose a choice of a clever basis for the representation of 
correlation structures for hierarchically organized stochastic processes: wavelets. They lead 



1 



to a tremendous simplification and directly unveil information of interest, which is otherwise 
scattered in the many dimensions of phase space. 



The call of wavelets ||10||-||14|| needs some further introductory remarks: The whole wavelet 
basis is constructed from dilations and translations of one single "mother" wavelet and rep- 
resents a selfsimilar and orthogonal basis. Thus we would expect that the wavelet transform 
uniquely simplifies the correlation functions of selfsimilar processes and, in particular, quasi- 
diagonalizes the covariance matrix. This has been demonstrated in ref. for fractional 



Brownian motion and in [|l(J for the p-model [||, which is the simplest tractable random 



cascade model. Still more can be learned from the wavelet transform. 

To estimate correlation patterns, the human brain follows a strategy different from the 
standard analysis of correlations: it organizes particles accumulated in densely populated 
regions into (hard to quantify) "clumps" or "clusters" and unpopulated regions into "voids" . 
If one looks closer into a particular "clump" it may (or may not) again be organized into 
"clusters" and "voids", but now with respect to the higher (smoothed) background density 
of the bigger "parent-clump", and so on. - Due to its multiscale localization property the 
wavelet transformation dissects a (random) signal into contributions from different scales 
and thus localizes small and large scale structures separately. This feature is the reason why 
wavelets are so successful in signal analysis [|TU . 



Wavelet correlations provide statistical information about further subclustering of certain 
structures like clusters, filaments or voids. Correlations between small subclusters living 
inside larger clusters or voids are naturally revealed in the wavelet transformed correlations, 
as are correlations between voids and filaments. Loosely speaking, we denote wavelet corre- 
lations as "clump" correlations. 

In this publication we study conventional and wavelet correlations of one-dimensional 
hierarchically organized branching processes. For those readers only interested in one or the 
other part, we have treated conventional and wavelet correlations in separated sections. 

In section 2 we define a class of discrete hierarchical branching models. We distinguish 
between the concepts of deterministic and random branching, which describe a zero lifetime 
cascade or a cascade with finite lifetime respectively. An evolution equation for the mul- 
tivariate generating function of the correlation densities is presented, which allows one to 
calculate the correlation densities recursively. This procedure is exemplified for the so-called 
p- model M, a-model H, the p- model with random branching and a QCD- motivated cascade 



model |T7|, [L8| . In an additional subsection we emphasize the conceptual distinction between 
an evolution (= theoretical) and a backward (= experimental) moment analysis. 

The concept of wavelet correlations is introduced in section 3. For the same models we 
discuss the simplified and compressed structure of the wavelet correlations of all orders. 
On the basis of the p-model we stress the relationship of higher order wavelet correlations 
with clump correlations. This interpretation is further illustrated as we discuss the wavelet 
multiresolution analysis of configurations belonging to two-dimensional branching processes. 

In section 4, which represents the conclusions, we give an outlook on future applications 
of the wavelet correlations in multiparticle processes and other fields in physics. 
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2 Hierarchical cascades 

— conventional correlation formalism 



We study a class of discrete hierarchical cascade models with what we call random or deter- 
ministic branching and restrict ourselves to the one-dimensional case only. From a general 
evolution equation for the generating function we calculate recursion relations for conven- 
tional correlation densities. As specific examples, the correlation structure of the frequently 
used p-model ||, ||] and a-model ||, [7j as well as the p-model with random branching and 
the QCD-motivated cascade model of ref. [|17], are analysed in more detail. We exhibit the 
importance of the conceptual distinction between moments referring to the forward evolution 
of the cascade and moments obtained from a backward analysis. 



2.1 Hierarchical cascading processes with random branching 

We illustrate the topological structure of a hierarchical cascading process with random 
branching in Fig.l. Starting from a single trunk, the tree might branch into a left and 
a right branch at the first cascade step with probability p; with probability 1 — p the tree 
might not branch at the first cascade step. At the second cascade step, every single branch 
that has formed so far may again split into a left and a right branch with probability p or 
remains unsplit with probability 1 — p. This prescription is repeated at each further cascade 
step. For lack of a better terminology we associate random branching with the case p < 1 
and term the case p = 1 deterministic branching. 

In addition to the topological structure of the branching process, the spatial structure is 
specified by the following rules (see Fig. 2): The origin of the cascading tree corresponds 
without any loss of generality to an interval [0, 1] with uniform (energy) density e normalized 
to one, = 1. If the interval does not split with probability 1 — p at the first cascade step, 
the energy density will remain the same as before. If, on the other hand, the interval does 
branch into two halves with probability p, one part of the original energy E — 1, namely 
E L = ^-E = goes to the left sub interval and another part, namely E R = goes to 
the right interval. The weights qi and q R are random variables in the range < q^, q^ < 2 
for both random and determinisitic branching and follow a joint probability distribution 
p(<1l,<1r), which is often called a "splitting function" or "splitting kernel". Depending on 
the splitting function p(qL, qn), energy may or may not be conserved at one branching step, 
so that El + Er = E needs not be fulfilled. The energy densities of the resulting left and 
right subintervals are = q^ and = q R respectively. - This prescription is repeated 
for all subsequent cascade steps j. Whenever a subinterval splits into two halves, the energy 
density of the left half subinterval is q^ times the previous energy density, e^k = QL^k^i 
and, correspondingly, the energy density of the right half subinterval is e^k+i^ = ( lR e k b ^ where 
qi and q^ always follow the same splitting function. The upper index of ej^ indicates the 
number of branchings jt, that occurred on the way to a subinterval and should not be confused 
with the number of cascade steps j, whereas the lower index numbers the subinterval given 
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at the branching scale If the subinterval does not split, the energy density remains the 
same as before. - This energy curdling is shown in Fig. 2 for one possible realization. Once 
the last cascade step j = J is reached, the resulting energy density distribution is resolved 
on this finest scale with 2 J subintervals (= bins). The energy densities e\f\ k=0,. . . ,2 J — 1, 
of neighbouring bins then may be identical, whenever a nonbranching has occurred during 
previous stages of the cascade. 

In particle physics, cascade models based on continuous random variables e k as discussed 
above are discretized by means of the Poisson transform mechanism ||. After the last 
cascade step J, the energy density ej^ is replaced by a discrete particle number n k which 
is drawn from a Poissonian with a mean proportional to ej^ /2 J . The necessary transition 
from ordinary to factorial correlation densities is well known, has no impact on our results 



and will not be further discussed in this paper. As has been demonstrated in ref. [|L6] this 
also holds for the transition from ordinary to factorial wavelet correlation densities. 



2.2 Evolution equation for conventional correlation densities 

The cascading process described in the previous subsection is selfsimilar by construction: At 
each cascade step the same branching prescription is applied as in the previous steps. As 
shown in later sections this selfsimilar construction does not necessarily imply perfect scale 
invariance of correlation functions. Nevertheless, a detailed analysis of correlation densities 
is the only means to obtain local information on the (selfsimilar) branching structure of a 
cascade process. 

The bin correlation densities are defined as the product of energy densities contained in 
various bins, averaged or sampled over all possible configurations. For example, one-, two- 
and three-bin correlation densities read 



(J) 
(J) 

Pk!k 2 

(J) 

Pkik-2k3 
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the brackets (...) indicate the averaging over all configurations. The indices fcj run from 
to 2 J — 1 and represent the number of different bins at the finest resolution scale J. 
The correlation densities are most easily determined once the corresponding (characterisitic) 
generating function 
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is known; they follow by taking appropriate derivatives with respect to the Ajj. . 
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In the following we develop a scheme how to determine the generating function at a given 
evolution scale from the previous evolution scales. This scale-dependent recursion relation 
then leads to scale-dependent recursion relations for the bin correlation densities. 

Suppose we know the generating function \v> after j cascade steps. Then we can 

determine Z^ +1 ^ after j + 1 cascade steps by a backward evolution: With probability (1 — 
p) there has been no branching in this backward evolution step and with probability p a 
branching into a left and right branch has occurred, where the split energy densities have 
been weighted with q L and q R respectively. This translates into 



Z U+D XU+D = (i -p) z® +p ldq L dq R p(q L ,q R ) 
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specify the spatial splitting structure of the branching process. The splitting function 
p(Ql, Qr) reflects the probability that random weights qL and q R are assigned to the energy 
densities of the left and right branch respectively and is normalized to one, i.e. / p(qL, Q R )dqLdq R 
1. The nonlinear backward evolution equation (f|) allows one to express Z^ +1 ^ in terms of 
Z^\ which in turn can be expressed in terms of Z^~ l \ and so on until the roughest scale 
j = is reached, where 
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The evolution equation (|J) has been derived in detail in |16| for the case of the p-model with 
deterministic branching. 

It is instructive to rewrite Eq.( 



as 
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where AZ^> 



A 



A 



Aj = 1 and 1/r = p/Aj is a branching rate. 



This form is a discrete analogue to the evolution equations widely used to model QCD 
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parton showers in the perturbative regime [19]. By this analogy we are motivated to study 
the mathematical structure of the recursive discrete evolution equation (|4|) in detail. 

From Eq. (|j) we derive the recursion relations for the bin correlation densities by taking 
the relevant derivatives (|3]) with respect to the A^ 1 ^ taking into account the transformations 
(H). For the one-bin correlation density we deduce: 

^ +1) = (i-p>K+p?pS ; (8) 

the following four cases have to be distinguished: 

• k\ G {0, . . . , 2 J ' — 1} = {L} and k\ = 2m an even number: 
m 1 = m = ki/2, m 2 = h,q = qZ = fp(q L ,qR) q L dq L dq R , 

• k± G {L} and k\ = 2m + 1 an odd number: 
mi = m = (ki - l)/2, m 2 = h, q = qZ, 

• ki G {2 J , . . . , 2 J+1 — 1} = {R} and k\ = 2m an even number: 
mi = m = ki/2, m 2 = h - 2\ q = q^ = J p(q L , Qr) Qr dq L dq R , 

• ki G {-R} and fci = 2m + 1 an odd number: 
m 1 = m = (k 1 - l)/2, m 2 = h- 2\ q = qk~. 

Because the splitting function p(g.L, g_R) needs not be symmetric in qi and q R we distinguish 
ql and ^r. For completeness we state that = 1, which follows directly from Eqs. (0) and 
- For the two-bin correlation density we have to distinguish two principal cases. If ki 
and k 2 belong to the same branch, i.e. ki, k 2 G {L} or ki, k 2 G {R}, we get 

pi + k ^(l-p)pS im2 +P^P ( S 3 m i (9) 

with q 2 = q\ = J p(qL, <?i?) q\ d^idg^ or q 2 = q 2 R analogously. If k\ and k 2 belong to different 
branches, i.e. ki G {L}, k 2 G {R} or vice versa, we get 

p£ 2 1} = (i - pjpffitn, + Pwmp^M (io) 

with gZ^R = / qn) qiqu^qL^qR- Several cases have to be distinguished for the indices 
ki-,k 2 — > m 1 ,m 2 ,m 3 ,m4, they follow from the transformations @ and will not be further 
specified. Again, for completeness, pjjj = 1. - Although the strategy how to calculate bin 
correlation densities of arbitrary order should be clear, we still want to mention three-bin 
correlations. Two principal cases have to be taken into account: if ki,k 2 , k% all belong to the 
same branch we have 

(i+i) = a_fi) n U) _L^^3 rt 0) n-n 

rkik2k:j V F)Fm\vri2mo, > F H rm^m^m% \ I 

with q 3 = q\ = J p(qi, qn) q\ &qiAqR or q 3 = q R , respectively, and if ki, k 2 belong for example 
to the left branch and k^ to the right branch we have 



(1 ~P)P + P QlQR Pm{m B Pm 6 ■ ( 12 ) 

The first value for the iteration is given by Pqoo = 1- 
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2.3 Specific cascade models 



The recursion relations @-(12) for the bin correlation densities are quite general. Models of 
hierarchical random branching processes differ in the branching probability and the splitting 
function, the latter giving rise to different splitting moments q^q 7 ^- In this subsection we 
concentrate first on two specific cascade models with determinstic branching, the p-model 
[||] and the a-model 0. They differ in that the former conserves energy at each branching 
whereas the latter does not. Both models have been suggested to describe the intermittent 
energy dissipation of fully developed turbulence; both also serve as simple discrete approxi- 
mations to multiparticle processes in high energy e + e~ and hadron-hadron collisions 0-0. 
Then we discuss the effects of random branching for the case of the p-model. Finally, a 
model in the context of perturbative QCD parton cascades [17], [18] is analyzed. 



2.3.1 p-model 

For the p-model |4| the branching probability is equal to p — 1; we call this deterministic 
branching. With every branching the energy splits into two unequal portions governed by 
the splitting factors = (1 + a)/2 and q R /2 = (1 — a)/2 or vice versa; both possibilities 
occur with probability 1/2. The splitting function of the p-model then reads: 

P(QL, q R ) = \(5{q L - (1 + a)) + 5(q L - (1 - a)))S(q L + q R - 2) . (13) 

The last 5-function in this expression guarantees that energy is conserved in the splitting 
process. Note that qi and q R are weight factors for the energy densities and not the energies, 
so that qi + q R = 2 and not equal to one. The moments of the splitting function which enter 
in the recursion relations for the correlation densities (§)-([l2|) now become: 

W = 1r = I p(Ql, qn) qi dq L dq R = 1 , 
Jo 

<ti = = 1 + « 2 > Wgn = f - a 2 , 

q\ = q 3 R = 1 + 3a 2 , q\q R = ~q~^ R = 1 - a 2 ; (14) 

they determine the scaling indices of the correlation densities. 

From the recursion relation @ it follows that the bin correlation densities of first order 
are equal to 

Pl? = l , (15) 

because the first order splitting moments in Eq. flT4D are equal to one. This is just another 
way of saying that on average the energy density at scale J is equal to the energy density at 
scale J = 0. For second order bin correlation densities we find from Eqs. (||), ([T0| ) and (|14D 

[ (i + «>il ,if h,k 2 e{L} , 

PkS = \ (i + ^U,^ ,Hh,k 2 e{R} , (16) 
^ (1 — a 2 ) , if k\ £ {L} and k 2 E {R} or vice versa. 
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The power-law-like singularity towards the diagonal can be seen clearly in Fig. 3a. The 
closer the bins are together topologically the more they share a common history and the 
stronger they are correlated. This is a characteristic feature of selfsimilarity. We shed 
some more light onto Eq. fll6|) , once we compare for example p$ = p^ = (1 + a 2 ) J and 
Pqi = Pio = (1 + a 2 ) J ~ 1 (l — a 2 ); they differ by a factor (l + a 2 )/(l — a 2 ). This represents an 
anticorrelation caused by energy conservation: 6q ^ and of one possible configuration 
differ insofar that only at the very last cascade step different energy-conserving weights 
(1 + a) and (1 — a) were assigned to them. 



2.3.2 a-model 



Like the p-model, the a- model is a deterministic cascade, so that p = 1. The splitting 
function factorizes and is now given as 

p(qL,qn) = (pi%L-(l-a))+P2%i-(l + /3))^ 

// ,-(l-a))+p 2 %H-(l + /3))) , (17) 



■(pi5(Qi 



where the splitting parameters a, (3 are positive real numbers. The probabilities p\ = 
(3 /{a + (3) and p 2 = a /(a + (3) are chosen such that p\ + p 2 = 1 and piil — a) +p 2 {l+ (3) = 1. 
In contrast to the p-model, energy needs not be conserved in every cascade step; for example 
the energy density may split into a left part with weight factor (1 — a) and into a right part 
with the same weight factor (1 — a) with a probability p\. Energy is however conserved on 
average since qT/R = 1. For the splitting moments appearing in eqs. fl8|)-(|T2D we derive from 



QE = Or = 1 , 

Ql = q 2 R = l + a(3 , qlM = 1 , 

~ql = 4 = l + af3(3-a + p) , ^ = ~q^q\ = 1 + a(3 . (18) 

With the first order splitting moments qZ, q^ and the branching probability p = 1 we 
derive from the recursion relation (^) 

P? = l • (19) 

On average the energy density contained in one bin is equal to the original energy density. 
According to the recursion relations @, (10) and the corresponding splitting moments 



( |T8"D we derive for the second order bin correlations 

= f (1 + a(3)p { £ k2 ,if kx, k 2 G {L} or . . . , . 

I klk2 \ 1 , if k x G {L} and k 2 e {R} or vice versa, 1 ; 

which are visualized in Fig. 3b. Apparently there is no big difference to the corresponding 
p-model relations ([TO]) : A powerlaw scaling towards the diagonal emerges, which is a clear 
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evidence of the underlying selfsimilarity in the a-model cascade. The difference to the 
p-model correlations lies in the absence of anticorrelations due to the violation of energy 
conservation. As soon as two different bins start to belong to different branches during the 
evolution of the cascade they become completely uncorrelated. We will see later on, that 
this missing anticorrelation is responsible for the deviation from perfect scaling behaviour of 
the a-model in a moment (multifractal) analysis. 

2.3.3 p-model with random branching 

In the previous subsections we have discussed the p- and a-model as two representatives 
of deterministic branching processes. Although these hierarchical models are intended as 
simple discrete approximations to multiparticle cascading @-[]7[|, they might not be that 
realistic. Multiparticle branchings have to be somehow "undeterministic" : With a certain 
probability a particle might not decay into two other particles at every step of the cascade; so 
to say the particle has a certain lifetime. Therefore we incorporate a nonbranching part into 
the hierarchical cascade process as described in sections 2.1 and 2.2. In fact, the branching 
rate 1/r introduced in the evolution equation (|7|) can be interpreted as the inverse of the 
particle's lifetime once the scale evolution parameter j is identified with time. With this 
motivation we study the effect of random branching on selfsimilar cascading. 

As in the deterministic p-model, we choose the splitting function ([Tj^) with the correspond- 
ing splitting moments ([14]). The branching probability p is now chosen in between < p < 1; 
the case p — 1 defines deterministic branching and the case p = describes no branching 
at all. The recursion relation (§) for the one-bin correlation density has as solution again 
p$ = 1. The solution of the recursion relations @ and (|10|) for the two-bin correlation 
density is depicted in Fig. 3c; there is an approximate power-law rise far from the diago- 
nal, but the closer the diagonal is approached the more the rise turns into a plateau. The 
explanation for this behaviour lies in the following: With an increasing number of cascade 
steps the probability that between neighbouring bins no branching has occurred becomes 
overwhelming. Then the energy densities of neighbouring bins are the same, so that the 
two-bin correlation density near the diagonal becomes constant. 



2.3.4 QCD-motivated cascade model 

As a generalization of the p- and a-model, so-called pseudo QCD-cascade models have been 
discussed in the context of multiparticle dynamics [PH| . Recently Brax, Meunier and Peschan- 
ski [|l7j] have suggested a discrete hierarchical branching model, which reproduces the same 
multiplicity moments as more elaborate gluon jet cascade models within perturbative QCD 
r^],[2l|-|23|. In this scenario the interval splitting at each branching step is interpreted as a 



partitioning in opening angle. The simplicity of the QCD motivated model of refs. []TT| , |TH| 
permits a full analytical treatment in terms of (local) correlation densities besides (global) 
moments. 
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The branching probability is set to p = 1, reflecting the "zero lifetime" of the virtual 
gluons. The splitting function reads 



p(Ql, Qr) 



5{q L - 2) S(q R ) + 





(21) 



and describes the decay of a hard parent gluon into one hard and one soft daughter gluon. 
The quantity (l/q) + = limp^o ((l/q)Q(q — 2(3) + S(q — 2(3) ln/3) indicates a regularization 

The parameter 7 

/aj. With (21) the splitting moments become 



prescription |2j 
strong coupling constant 



7oln2 = w6a s /7rln2 is proportional to the fixed 



Ql 



Qr 



Qr 



Qr 



dq L dqRp(qL,q R )qL =1 + 7 

qEqr = 27 



2 1 + 



7^ 



4(1 + 



QlQr 



QlQr = 67 



(22) 



This sets the stage for the conventional correlation densities (jq)- (|12l) 
For the first order correlation density we find 



(J) 
Pki 



;i+7)* 



(23) 



The splitting function ([H]) does not conserve energy; it produces additional energy in every 
cascade step. Here, we should think of "energy" as "multiplicity", so that ( p3|) describes the 
average production of multiplicity during the evolution of the cascade. 

With the help of the recursion relations @ and fllOD , the relevant splitting moments ( p2|) 
and the result (|23|) we calculate for the bin correlation densities of second order 



.Cj+1) 

Pk x k 2 




l)pk!k 2 
2i 



, if ki, hi G {L} or . . . 

, if ki € {L} and k2 G {R} or vice versa. 



(24) 



The reduced second order correlation density pu^l PkiPk 2 * s shown in Fig. 3d. Note that for 
the p- and a-model with and without random branching it was superfluous to distinguish 
between a reduced and non-reduced correlation density because there pu = 1. This is not 
anymore the case for the QCD-motivated cascade model treated here. 

From Eq. (g4|) it follows with fci = that p Q ±/po,2k = (2 + 7)/(l + 7) 2 for k ^ and 
Po,o/po,i = (2 + 7V27. This implies a scaling [(2 + 7)/(l +7) 2 p towards the diagonal except 
for the last step. It is the same scaling we encounter in a (forward) normalized moment 
analysis in section 2.4.2. Due to the pronounced "non-scaling" peak along the diagonal and 
the non-matching anticorrelations, the correct (backward) moments will deviate from this 
perfect scaling. 
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2.4 Moments: forward evolution vs. backward analysis 

Contrary to common belief, the a-model does not exhibit a perfect power law in a (backward) 
moment analysis in spite of its strictly selfsimilar (forward) evolution. This statement also 
holds for the QCD-motivated cascade model discussed above. In order to unveil this hidden 
secret we recall the standard procedure of a moment analysis as proposed in ref. f|. 

For example, the second order moment M 2 (J, J) is defined by 

m 2 ( j, j) = i 2 £ ({A J) f) = ^ 2 £ p { 2 = p\£ , (25) 

where the first index j x in the nonstandard, double- indexed notation M 2 (j%, j 2 ) denotes the 
total number of cascade steps J after a forward evolution whereas the second index j 2 < j\ 
represents the resolution scale (~ log(bin size)), at which the cascade is analyzed in backward 
direction. Note that the bin correlation densities pj^. of the p- and a-model as well as of the 
QCD-motivated cascade model are independent of the shift index k. Now we average over 
two neighbouring bins and arrive at the moment 

m 2 (j, j - 1) = ^ 2 e 1 ( Q(4£ + O) 2 ) = \{p { £ + P ( oi ] ) = 1(m 2 (j, j) + p^) , 

(26) 

resolved at the scale j 2 = J—l. Then, we continue to average over four, eight, . . . neighbouring 
bins and arrive at the following recursion relation for the moments: 

M 2 (J, J - j - 1) = ~ (M 2 ( J, J — j) + p { %) . (27) 
2.4.1 p- and a-model 

For the p-model we realize with the help of ( |i~6D that 

p-model: M 2 (J, j) = M 2 (j, j) = p® = (1 + a 2 y . (28) 

The moment obtained after J cascade steps, averaged over bins of size 2 J ~- 7 and then again 
averaged over 2 J blocks of bins is equal to the diagonal element of the bin correlation density 
obtained after j cascade steps. 

This does not hold anymore for the a-model: From (^7|) and (|20|) we derive 

a-model: M 2 (J,j) = - 1 - (1 + a(3) j \ 1 - af3 

(1 - a[3) y 

In Fig. 4 we summarize the various scaling aspects of the moments M 2 (J,j). Obviously, 
the backward analyzed moment M 2 (J,j) is not equal to the bin correlation density p^l = 



(^n • <29) 
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(1 + a/3y = M 2 (j,j) obtained after j forward a-model cascade steps. The reason for this lies 
in the fact that energy is not conserved in the a-model; an energy density ejjf obtained after j 
cascade steps might split into two energy densities e 2 ^ = (l+/3)ej^ and eik+l = (l+/^) e A? ^ m 
the next cascade step. But then ejjf 5 ^ §(^ +1) + 4fc+?) so tliat M 2 (j + l,j) ^ pjjjj = M 2 (j,j). 

Here it is really important to distinguish between the moment analysis in evolution direc- 
tion of the cascade, i.e. m 2 {j) = M 2 (j,j) in dependence of j, and in backward direction after 
J cascade steps, i.e. M 2 (J,j) in dependence of j. An experimental moment analysis is always 
based on moments like M 2 (J,j): every configuration (event) is analyzed backwards, given a 
fully developed cascade at final scale J. From these considerations about the a-model we 
learn, that the (observed) moments M 2 (J,j) may not be naively identified with the forward 
evolution moments m 2 (j). This is in contrast with the results of the p-model cascade, where 
the anticorrelations due to energy conservation make a distinction between an evolution and 
a backward analysis unnecessary. 



2.4.2 QCD-motivated cascade model 



From the (backward) moment recursion relation (|2~7 ) and Eq. ( p4]) we derive for the QCD- 
motivated cascade model 

M 2 (j,j) = 2 (1 + ^ (*±r) S + l±*L (i + 1) J v . (30) 

V ,JJ 3 + 2 7 V(1 + 7)V 3 + 2 7 V 2J y ! 

This is an unnormalized moment. For a proper normalization we have to divide by the 
square of the density (p[f^) 2 = (1 + l) 2J from Eq. (^); this normalization does not affect 
the scale j dependence. 

We observe that the above expression consists of two terms with different scaling behaviour. 
The first term scales with a factor we have already noticed in the second order correlation 
densities (p4|), whereas the second term scales trivially in the factor 2. This can be explained 
in the following way: The second order correlation density ( p4|) can be decomposed into a 
correlation matrix with exact anticorrelations and an additional purely diagonal correlation 
matrix; the former gives rise to the first term in (|3l]) with anomalous scaling, whereas the 
latter accounts for the trivial scaling in powers of 2 of the second term. 

Note in this context, that the expression for the second order moment (]29f) of the a-model 
can also be cast in this form: one term scales with the factor (1 + aft), whereas the second 
term scales with the factor 2. In contrast to the a-model result, the second term in ( p0|) 
for the QCD-motivated cascade model comes with a positive sign. As a consequence, in a 
log-log plot the curve for M 2 bends upwards for the fine resolution scales, i.e. large j values, 
and deviates from the perfect scaling shown at the rougher resolution scales, i.e. small j 
values, where the first term in (^) dominates over the second one. This result is illustrated 
in Fig. 5; compare also again with the corresponding a-model result shown in Fig. 4. 

We emphasize again that the (backward) moment M 2 (J,j) of (|30| ) is the relevant moment 
to be compared to an experimental analysis. In theory, the normalized moments m q (j) 
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in forward evolution of the cascade are calculated in most cases and compared directly to 
experimental data. For example, the corresponding normalized moment of second order is 
defined as 

Here the diagonal elements of the correlation densities enter for a cascade with j cascade 
steps; the cascade is not evolved down to the finest scale J and then resolved backwards on 
a rougher resolution scale j. With the results (^) and (|23j) we find for the QCD- motivated 
cascade model the scaling law 



m 2 (j) 



2 



2 + 7 V 
(1+7)V 

log 2 (2+ 7 )+21og 2 (l+ 7 ) 



-J 



V^r (1 ~ |7o) (for 7o «l). (32) 



It shows the same anomalous scaling as the first term in the corresponding (backward) 
moment ( |30] ) and is also illustrated in Fig. 5. The exponent |7 appearing in the last 
expression is identified with the Renyi dimension D 2 . 

For the (forward) normalized moments of higher order the calculation is straightforward; 
we state the outcome without giving further details: 



m q 



^ _ ^ 2 _^-log 2 [2 9 - 1 (l+^)]+glog 2 (l+ 7 ) 

(for 7o <l) (33) 



with 

= — 7b . (34) 

q 

These scaling exponents for the (forward) evolution moments have been derived previously 
in the context of a perturbative QCD gluon cascade [ ITU , 1 21 ]— [ 231 - Note again that it is not 



these unobservable (forward) evolution moments which should be compared to experimental 
data, but only experimentally accessible (backward) moments. 



3 Wavelet correlations in hierarchical branching pro- 
cesses 

The conventional correlation densities (HD characterize the correlations between bin contents 
resolved at some finest scale J (usually chosen somewhere near the resolution limit of the 
measurement device). They do not provide a natural description in terms of larger structures 
(clusters and voids). However, sometimes a choice of a clever basis for the representation 
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of the bin contents can directly unveil this information of interest. The wavelet transform 
10|| -[ l4ll appears to be such a clever change of basis functions for hierarchically organized 



stochastic processes JT3], The strategy applied in this section is the following: We 

apply a multiresolution analysis arising from the wavelet transform to each realization of 
the branching process and study the correlations between the wavelet amplitudes. This 
procedure is exemplified for the same cascade models discussed in section 2. 



3.1 Multiresolution analysis 

In this subsection we will briefly recite only the most important ingredients of a wavelet- 
based multiresolution analysis. We will concentrate on the Haar wavelet only, which is the 
simplest of all wavelets; generalizations to other wavelets are straightforward. 

The binned energy densities e£ can be viewed as a step function of x G [0, 1]: 



where 



e (J) (x) = £ ei J Ujk(x) , (35) 

fe=0 



a < \ ^<oJ 1A J 1 for k2 ~ J <x<(k + 1)2- J , 
<pjk{x) = 0(2 x - k) = i Q ^ (36) 

are simple box functions. In this respect the energy densities e k J ^ can be interpreted as 
amplitudes of the orthogonal expansion of e^ J ^(x) in terms of the basis functions 4>,jk{x) at 
scale J. 

The fundamental dilation equations 

<t>j-\,k{x) = 0(2 J ~ 1 x - k) = £ c m <f){^x -2k-m) 

m 

— ^^C m (j)j t 2k+m{.x) , 
m 

^-i, k {x) = ^{2 j - l x-k) = J2(-~L) rn ci- m <P{2 j x-2k-m) 

m 

= ^(-l) m Cl _ m Ji2fc+m (x) (37) 

m 

define a multiresolution analysis |T0|-[r2|. The associated functions <ft and ip are called scaling 
function and wavelet, respectively. For the box functions (p6|) we have Cq — c\ — 1 and = 
otherwise. A different choice for admissible coefficients Ck corresponds to a different filter for 
the multiresolution analysis and leads to a different wavelet ip{x)] see ref. |yj for a variety 
of examples. 

In Eq. fl3"?D the box functions (f)jk(x) at scale J are expressed in terms of the box functions 
(fij-i t k'(x) at scale J — 1 and the "difference functions" iftj-i,k"(x) at scale J — 1 with 
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The set of functions ipjk{%) is called Haar- wavelet basis. Further iterating Eqs. fl3~T|), the 
energy density step function J ) (x) of Eq. ( j35"l) can be expressed as a double sum 

e (J) (x) = EE^'^W + 4°Voo(x) (39) 

3=0 fc=0 

representing the multiscale decomposition. In other words, e^ J '{x) is dissected into contri- 
butions from different scales j. Eq. (^) defines the linear wavelet transformation, which is 
governed by the coefficients c k entering Eq. (|3T|): 



with 



|=We = W(c fc )e (40) 



e = [£o, ei, • • • , £2J-i) , 

- , (0) ~ ~ s 

e — l e o > e oo, eio, en, £20, • • • , 

=: (e , ei, e 2 , . . . , e 2 j_i) , (41) 

where the upper index (J) has been omitted. The transformed amplitudes e*-^ are called 
wavelet amplitudes. The explicit form of the transformation matrix W is exhibited for 
example in ref. [Of . 



3.2 Wavelet transformed correlation densities 

Correlations between the transformed amplitudes are called wavelet correlations and have 
been introduced already in previous publications |I5|, [HJ. The wavelet correlation densities 
follow by applying the wavelet transformation either directly to the conventional correlation 
densities or to the evolution equation for the corresponding generating function. 

If the bin correlations ([J) are known, the wavelet correlations can be directly deduced with 
the help of ©: 

pfcx = (hi) = ( W fclfe2 e fc2 \ = ^ W fclfc2 (e k2 ) = ^2W klk2 p k2 , 

Pkikz = (hih 3 ) = ^2 ^kik2^ r k 3 k 4 Pk 2 k 4 , 

k 2 k 4 

Pkik 3 k 5 = (hih 3 h B ) = X! Wfcj^W^^W^^p^fc^u . (42) 

k 2 k 4 k 6 

Thus the wavelet correlations can be obtained directly from the standard bin correlations. 
The treatment of wavelet correlations is general and does not depend on the specific choice 
of admissible coefficients c k defining different wavelets. 

For the case of the simple Haar wavelet it is instructive to deduce the wavelet correlations 
directly from the generating function (§). This gives more analytical insight into the corre- 
lation structure especially for the binary branching processes we have introduced in section 
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2. In the exponent of the generating function (|2]) we introduce the change into the Haar 
wavelet basis such that 



2 J -1 

E A fcl e fcl 

fci=0 



2 J -1 2 J -1 

fci,/C3=0 fel,fc2,fc3=0 
2 J -1 2 J -1 

E (W _1 ) fcl fc 2 ^fciefc 2 = E %2^fc 2 • 

ki,k2=0 &2=0 



(43) 



The coordinates r/k show a simple relation with respect to the A^: 



2 J -1 



E ^ > 

k=0 







2 J ~ 1 -1 


2 J -1 






m 


= Voo 


= E A fe - 


E a, 


1 








k=0 


fc=2 J -! 










2 J ~ 2 -l 




3 . 2 J-2_ 1 


2 J -1 


V2 


= Vw 


= E a,- 


E a, 


= VU = E X k ~ 


E A fc 






k=0 


fc=2 J ~ 2 


k=2 J - 1 


fc=3-2 J - 2 


7] 2 J-i 


= Vj- 


1,0 — Aq — Ai 


, . . . , 


V2 J -1 = Vj-l^J- 1 -! = A 2 J-2 " 


- -V-i ; 



(44) 



for more details see again ref . JTB] . In view of the evolution equation ([|) with the splitting 
rules (H) it is necessary to express ff-^ +r> at scale j+1 in terms of ffi\ ff^ and fflj at scale 
j; we deduce: 



„(i+i) Aj+i) _C?+i) _0"+i) fj+i) „C?+i) 

70 ) 700 5 710 ) '/ll 5 720 ) •••) 'lj,2i-lj ' 

,0+i) , „0+i) „0+i) Jj+i) „0'+i) 



■l(V) U+ij 4- f ,^ +1 ^ r) U+ij r) U+ij r ) U + 1 ' 1 n U+iJ r) (i+1) 
2^/0 "T 700 J; 710 > '120 j '121 j 730 ) •••! 'IjflJ- 1 -! 

( 0) 0) 0) 0) 0) 0) 

l//0;L' '/00;L' 710;£> 711;£j 720;L> • • • ' 'A/-l l2 i-l-l 



LI ' 



-n U+1) \ n ij+1) n {j+1) n U+1) 

2"v'0 700 J' 711 5 722 ) 723 ) 734 

f Jj) JJ) JJ) Jj) Jj) Jj) , 

(//0;.R> 700;ii5 710;R) 711;iiJ 720;i?> •••> 'lj-l,2i- l -l;Rj ' 

(J+1) 
' 7j-l )2 J-l_l 

Si) 



i ■ ■ ■ i Vj^j-i 



0+i) 



_CM-i) _(?+!) 

\//0 > 700 5 710 ) 711 ) 720 5 

( 0) 0) 0) 0) 0) 

\\0\Mi 700;M) 710;M' 711;M' 720;AfJ • • • ) '\j-\,2i- x -\;M 



(45) 



This sets the stage for the evolution equation of the generating function for the Haar wavelet 
correlations; from Eq. M) we deduce straightforwardly 



^•(i+i) 



^O+i) 



1 -p)Z^[#] +p f dq L dq R p(q L ,q R ) Z^[q L rfi ] ] 



QrVr 



(46) 
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Taking the derivatives with respect to the Vjiti \ recursion relations for the Haar wavelet 
correlations can be found. For the first order Haar wavelet correlations we find 



Pt 1] 



5 U+1) 

P(00) 
3CJ+I) 



1 dZU+V[ffj+ 1 )} 



dm 



(j+i) 



»?(i+ 1 )=0 



:i-p)(i-Mptoi)+pV^# 



(47) 



the last relation holds for 1 < ji < j and q = qi, & 2 = k\ for (ji^i) G {L} or 5 = g^, 
k 2 = ki — 2 J1_1 for (jifci) G {-R}. The starting value is = p^ = 1. Furthermore, we 
specify three recursion relations for the second order Haar wavelet correlations, which will 
turn out to be important for later discussions: 



Po,o 



i_ a 2 zk' +i y +i >] 

a \V0 ) ff 3 + 1) =0 

( 1 - P)Po,o + P j Po,o + P 



4 



~(i+l) 

"(00), (00) 

^Oifci),(iifci) 



where again g 2 = q\ or and /c 2 = k\ or fci 
relation. 



(! -P)(! - ^o)P(oo),(oo) +P -^-Pofl-P—iPo ) ' 

(1 -p)(l - ^•JpgJ^.^fa) +P^pS-i,fo),(ii-i,fe) > ( 48 ) 
_ 2 jW for (jifci) G {L} or {#} in the last 



3.3 Wavelet correlations of specific cascade models 
3.3.1 p-model 



Wavelet correlations in the p-model were already studied in ref. fig] , emphasizing the self- 
similarity aspect of the wavelet basis. We briefly recall the main results. Since the splitting 



moments q^ and g# are equal (see Eq. fll4|) ), the first order Haar wavelet correlations ( f47|) 
are all zero, 

/5(j*) = , (49) 

except for po = 1. On average a difference between energy densities of neighbouring bins 
is zero, as a difference might come with a positive or negative sign with equal probability. 
Only the global average po of all bins is one. 

From (HST) and (O) we deduce for the second order Haar wavelet correlations 



Po,o = 1 , 

P(ji*i).(ia*a) = a2 ( X + a2 y is hh s kik 2 ; (50) 
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they are depicted in Fig. 6a, where the wavelet indices (jk) are ordered according to Eq. ([41]). 
The product of two differences belonging to different scales or bin positions is zero on average. 
The diagonal contributions, which represent squares of differences, show a scaling law as the 



resolution ji increases; compare with (|28|) . It is exactly this diagonal structure we would have 
expected from second order Haar wavelet correlations of the selfsimilar p-model cascade. In 
ref. |16| we have shown that when choosing wavelets other than the Haar wavelet, the second 
order wavelet correlations of the p-model cascade turn out to be quasidiagonal, with diagonal 
contributions dominating over the off-diagonal ones and also showing an approximate scaling 
law. 

For the Haar wavelet correlations of the third order a double scaling is found: 

/SoifcO.CfcfaO.Cfrfe.) = (1 + 3a 2 ) J1 (±2« 4 (1 + a 2 )'" 3 "*- 1 ) , (51) 

where (jik\), C/2^2) have to share a common parenthood, i.e. < j\ < A; 1 2 jf2-Jl < k 2 < 
k{2?*-n + - 1 for the (+) sign and k{2?^ + 2^ 1 ~ 1 < k 2 < (k x + 1)2^-^ - 1 for 

the (-) sign. All the other pure third order Haar wavelet correlations are zero; only some of 
those involving an index 0, which represents the amplitude e = 1 of the global average and 
thus not a difference amplitude, are nonzero but would vanish for Haar wavelet cumulants 

g. 



Note that the wavelet transform compresses the full information contained in the second 
order correlation function ( |1~6"D into the diagonal fl5"0|); compare also Fig. 4a with Fig. 6a. 
In addition to the diagonal contributions also certain off-diagonal bandstructures arise in 
higher order wavelet correlations as e.g. in Eq. (^TJ). An interpretation of these offdiagonal 
contributions in terms of clump correlations is given in section 3.4. 



3.3.2 a-model 



From Eqs. (|47|) and (|T8|) we derive exactly the same result for the first order Haar wavelet 
correlations of the a-model as stated in eq. (^) for the p-model. Differences between the 
two arise for higher order Haar wavelet correlations. With the branching probability p = 1 
and the second order splitting moments of (fHJ), we derive 



~(j) 
Po,o 



1 



'l-a(3) 




_ (J) _ qg(l + a/?)* ) ( l + aP \ J j \ , 

P(hki),(j2k 2 ) ~ ( X _ a £) 1 1 I 9 ] f °hh°kik 2 [ZA) 

from the recursion relations (|4"8| ) for the second order Haar wavelet correlations. This corre- 
lation matrix, which is depicted in Fig. 6b, is diagonal as in the p-model case. However, the 
diagonal contributions now depend on the number of performed cascade steps J. 

We first comment on the element pgci- m contrast to the p- model result, it is unequal to 
one; for \a(3\ < 1 it is equal to 1/(1 — a(3) in the limit J — > 00. The quantity poci — 1 reflects 
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the square width of the fluctuation in the total energy. Its deviation from zero is a clear 
evidence that the total energy is not conserved in the a-model cascade. 

The elements P(hki),(jiki)> reflecting the true second order Haar wavelet correlations, do 
not show a rigorous scaling; only for ji values being "safely" smaller than the number of 
cascade steps J and |a/3| < 1 does a power law depending on the scale index ji show up. If 
we set a = (3 in the a-model and use the same value as in the p-model, we obtain the same 
asymptotic scaling as for the perfect scaling in the p-model. The only difference is that the 
missing anticorrelations lead to a modification of 1/(1 — a(3) in the overall factor. 

In Fig. 4, we show the deviation from scaling of the diagonal elements of the second order 
Haar wavelet transformed correlations p/\ ^ with respect to the scale index j for both 
the a- and the p-model. Note the resemblance between Eq. fl52[) and the moment Eq. (p9f), 
which emphasizes again that the wavelet correlations imply, by design, a backward analysis. 
Note also by comparison of Figs. 3b and 6b that these deviations from perfect scaling are 
observed more easily in the wavelet picture. 



3.3.3 p-model with random branching 



We do not learn anything new from the recursion relations (17) for the first order Haar 



wavelet correlations. The global average is p^ = 1; first order differences vanish for any p, 
so that pijl} = 0. 



Turning to the more interesting second order Haar wavelet correlations (]48|) , we find, as 
expected, that Pqq = 1. Again, energy is conserved in every possible configuration. Off- 
diagonal contributions to the second order Haar wavelet correlation matrix pS^) fefo) vamsn 
also in this more general cascade. Next, we consider the two diagonal elements p^ Q L ^ ^ and 
in particular. 

For P(o,o),(o,o) we find 

P(o,o),(o,o)=« 2 (l-(l-P) J ) • (53) 

The factor (1 — p) J represents the probability that after J cascade steps no branching has 
occurred. For < p < 1 and J — > oo this probability goes to zero. Then ff/p \ / y which 
represents the square of the difference between the average of the energy densities belonging 
to the 2 J_1 left bins and the average of those belonging to the 2 J_1 right bins, approaches 
the value (|50|) obtained for the p- model with p — 1. 

The element pfj^ k \ k \ representing the square of the difference between energy densi- 
ties in two adjacent bins at the finest scale follows from the third equation of fl4~8|): 

j-i 



P(j-i,kW-i,k)=P* 2 (p(l + * 2 )) ~ • (54) 



For p < 1/(1 + a 2 ) this element tends to zero as J —>■ oo. As the number of cascade steps 
increases, the probability that no branching has occurred on the finest scales also increases, so 
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that the difference amplitudes on the finest scales will be zero. Once p > 1/ (1 + a 2 ) however, 

this interpretation no longer holds. For p = 1/(1 + a 2 ) we have always pfj-i^),u-i,k) = P a2 i 
no matter how large the number of cascade steps J. Somehow this seems to contradict 
common sense, especially when the number of cascade steps becomes very large. It is true 
that, as J — > oo, the probability of a branching occurring at every cascade step becomes 
infinitesimal small; but on the other hand, once it happens, the difference amplitudes on 
the finest scales become very large because neighbouring bins have gone through a long 
common history and have built up huge energy densities. Differences of large densities are 
also large since the splitting in the very last cascade step always goes with a (1 + a) to 
the left and a (1 — a) to the right or vice versa. This huge difference then dominates over 
the small probability that branching occurs; as a consequence the correlation p^j_ x k \ rj_ x k \ 
stays finite and unequal zero for p = 1/(1 + a 2 ). For p > 1/(1 + a 2 ) it even increases as the 
number of cascade steps J goes to infinity. 

Fig. 6c shows the full second order Haar wavelet correlation density for the p-model with 
random branching for the case p < 1/(1 + a 2 ). In Fig. 7, the scale dependence of Pu\\uk) ^ s 
shown for various values of the branching probability p. For j much smaller than J the same 
scaling shows up for the p-model with random branching as with deterministic branching; 
large-scale branchings have occurred during the cascade with a probability almost equal to 
one. For small branching probabilities, i.e. p < 1/(1 + a 2 ), a deviation from scaling sets 
in at midscales and the contributions p$ k \ (j k ) drop rapidly to zero as finer and finer scales 
are considered. This is clearly a lifetime effect, because branchings at the finest scales will 
hardly occur. 

For higher order Haar wavelet correlations similar modifications are to be expected. The 
multiple scalings ( |5TD and (|60| ) for higher order Haar wavelet correlations will only live on 
the larger scales and will break down at the finest scales. 



3.3.4 QCD-motivated cascade model 

The first order Haar wavelet correlation densities for the QCD-motivated cascade model 
follow from the recursion relation ( fjT|) using the splitting moments (p2|) and p = 1: 

A J) = (1+7) J , 

As expected, the global average multiplicity density p^ is equal to the local average density 
p k of Eq. (p3|) , whereas all local differences are zero on average. 



From ( flq) and (p22|) we deduce for the second order Haar wavelet correlations 

-M _ (1+7)" , (, , 7- J 
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~(J) 



(1-7-7 2 ) 



+ 



1+7 
(1 + 7 
(1 + ^) 



7) 



2J-2 



1 + 



7 



2 J 



The average square of the total multiplicity density pofo i s unequal to {po 1 ') 2 , signaling 
the absence of proper anticorrelations. The true second order Haar wavelet correlations 
PUiki) (j2k 2 ) are a S am diagonal; normalized with respect to (pb^) 2 they are shown in Fig. 6d. 
The diagonal contributions depend on the number of performed cascade steps J and do not 
show rigorous scaling. In fact, as in the case of (backward) moments (p0|), the dependence 
on the resolution scale j splits into a term with the same anomalous scaling and another 
term scaling trivially in powers of 2. This j-dependence is once more illustrated in Fig. 5; 
the deviation from the anomalous scaling at the rougher resolution scales sets in earlier for 
the wavelet correlations than for the (backward) moments. 

It is instructive to introduce wavelet moments W q (J,j). For the case of second order we 
define 

^(J,j) = ljEP { uku k) = $U) • (57) 

z fc=0 

They are related to the (backward) moments M2(J,j) of ( E7J ) via 



2 + 7 
(l + 7) ; 



(56) 



W 2 {J,J-j-l) 



M 2 (J,J-j)-M 2 (J,J-j 
l -(M 2 {J,J-j) 



(J) 



(5f 



which follows from Eq. (42|). This relation reflects the multiresolution property of the 
wavelets as they look on the difference of two adjacent scales. A "difference" moment is 
more sensitive to deviations from perfect scaling than an "average" moment; this explains 
the results depicted in Fig. 5. Also, in view of the relation fl58|), it now becomes clear that 
the expression (|56"D for the wavelet correlations splits into the same two different scaling 
terms as the (backward) moments fl3"0|) did. 



Expressions for higher order wavelet correlation densities will not be given here for the 
QCD-motivated cascade model. Their general interpretation will be given in the next two 
sections 3.4 and 3.5. For a gluon cascade they can be understood as correlations between 
gluon subjets inside larger gluon jets. 



3.4 Higher order wavelet correlations: clumps 

What does the double scaling of the third order wavelet correlation fl5lD tell us? To get 
a glimpse, we consider once again the evolution of the energy densities according to the 
cascading prescription of the |?-model. See Fig. 8. After j 1 cascade steps we pick out the 
energy density ej£ of one bin with label k\. At the next cascade step this bin has split 
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into two subbins with labels 2ki and 2ki + 1 respectively, which contain the energy densities 
e 2fci +1 ' ) an d 4irfi ■ We ask now for the correlation between the energy density ejf 2 contained 
in a subbin with labels j 2 , k 2 and the energy density ej^ contained in the picked bin with 
labels ji,ki where the support of the subbin j 2 , k 2 is contained in the support of the bin 
ji, k\\ i.e. j\ < j 2 and fci2 3 ' 2 "" 3 ' 1 < k 2 < {k\ + l)2 jf2_J1 — 1. We get for example: 

(^ 1) 4 2) > = ((^ ) ) 2 ><*" J1) > = (1 + « 2 ) J1 . 
^ 1 1) (4 2) ) 2 > = ((4 l) ) 3 ><(4 2 ~ Jl) ) 2 ) = ( 1 + 3« 2 ) J1 (l+« 2 )^ 1 , (59) 

where k 3 = k 2 — /ci2 J,2_J1 . These interscale correlations contain additional information about 
the subclustering structure of the cascade process. In other words, this can be interpreted 
as the correlation between a "coherent structure" and one of its "substructures". Such 
(sub) structures can be identified with clusters (high density regions) as well as voids (low 
density regions) and will subsequently be referred to as "clumps" . 

The second order correlation (ej^ej^ ) only provides information about the common 
"history" of the clumps j\kx and j 2 k 2 and not about the subclustering structure, because 
( e k3 = 1- This situation changes once we consider the third order correlation (ejf 1 1 ' ) (ejf 2 2 ' ) ) 2 ). 
Here we get a double scaling, which depends on the common "history" of the two clumps as 
well as on the subclustering structure of the larger clump. It is the same double scaling we 
have found for the third order Haar wavelet correlations ( |5"TD . Therefore we interpret higher 
order Haar wavelet correlations as clump correlations. In fact, what the third order Haar 
wavelet correlations do, is that they correlate the difference of the energy densities of two 
adjacent clumps living on the scale j\ with the square of the difference of the energy densities 
of two neighbouring sub clumps living on the scale j 2 inside the original clump, as depicted 
in Fig. 8 by marked arrows. Of course, this is just another definition of clump correlations. 

For completeness we state the nonvanishing contributions to the fourth order Haar wavelet 
correlations of the p-model which provide further information about the subclumping struc- 
ture: 

P(iifci) 4 = a4 (l + 6 « 2 + a 4 ) J1 , 
p (ilfcl)2)(i2fc2)2 = a\l + a 2 r^{l + Qa 2 + a A r , 

(j2 < ji, same parenthood), 
= a\l - a 2 )\l + 6a 2 + a 4 ) s (l + a 2 )^' 28 ' 2 , 

(partially same parenthood), 

PWWiUisks) = (±) h (±) h 2« 6 (3 + a 2 )(l + 6a 2 + a 4 )^(l + 3a 3 )*-*- 1 • 

.(1 + a yi-h-i ; 

(j3 < 32 < jii same parenthood). (60) 

The phrase "same parenthood" used in the first case of Pu 1 k 1 ) 2 ,(j2k 2 ) 2 stands for Z^- 71- - 72 < 
ki < {k 2 + 1)2 J ' 1-J ' 2 — 1. The second case with "partially same parenthood" translates into 
< k < 2 s - 1, k ■ 2^- s <h<(k + \)2^- s - 1, (k + |)2^- s < k 2 < (k + l)2^~ s - 1, where 
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s represents the scale, from which point the wavelet indices (ji&a) and {j 2 k 2 ) follow different 
branches in the underlying tree structure, (ji — s) and (j 2 — s) are the scales of two different 
subclumps relative to the scale s of their common parent clump. Thus the correlation density 
P(j 1 k 1 ) 2 ,{j2k 2 ) 2 provides information not only about the direct subclustering of a large clump, 
but also relates two different subclumps within a large clump. 

Furthermore, the triple scaling of P(j 1 ki) 2 ,(j2k2),{j3k 3 ) fells us something about the correlation 
of a subsubclump within a subclump within a clump. Here the phrase "same parenthood" 
illustrates k 3 ■ 2 j2 ~ h < k 2 < (k 3 + \)Z j2 ~ j3 - 1 for the (+) sign and (k 3 + \)2 j2 ~^ < k 2 < 
(fc 3 + l)2^ 3 -l for the (-) sign in (±) h and, furthermore, k 2 -2^^ 2 <k x < (k 2 + \)2^i 2 -I 
for the (+) sign and (k 2 + |)2^ 2 < k x < (k 2 + 1)2^ 2 - 1 for the (-) sign in (±) h . - This 
demonstrates once more that higher order wavelet correlations provide direct information 
about the clustering hierarchy and can be interpreted as clump correlations. In the next 
section we further illuminate this interpretation from a two-dimensional perspective. 

We emphasize that this kind of subclustering information cannot be obtained by a conven- 
tional moment analysis along the lines of section 2.4. Although in principle this information 
is contained in the conventional higher-order correlation densities, it is very hard to measure 
or visualize directly. Due to its multiresolution character, the wavelet transform compresses 
exactly this information of interest to a readily accessible form. 

3.5 Multiscale clustering 

in two-dimensional branching processes 

To illustrate and visualize clumps in more detail, we will now consider hierarchical branching 
processes in two dimensions. As a representative we choose the two-dimensional a-model. 
Compared to the p-model, the a-model is easier to generalize to higher dimensions, as it can 
be built by direct products of one- dimensional a-models. Instead of halving a given interval, 
a square is subdivided into four subsquares with labels 1,. . . ,4 during one cascade step. The 
original energy is partitioned according to the splitting function 

p(?i, , 94) = II - (i - «)) + pA<h - (1 + /?))] ; (6i) 

i=l 

this is a straightforward generalisation of eq. (|T7[) . The prescription ( |61[ ) is repeated for all 
following cascade steps. One possible realization of the two-dimensional a-model is depicted 
in the upper left corner of Fig. 9; six cascade steps (j = 6) have been performed with the 
splitting parameters a — (3 — 0.4 We have used a gray scale to indicate the population of 
regions in between large energy densities (white) and small energy densities (black). We 
observe that certain regions clump into clusters of large densities and other regions with 
small energy densities exist as voids. Furthermore clusters/ voids appear in different sizes 
and show substructures, which are again clumped into clusters and voids. 

To quantify this picture, we explicitly perform a multiresolution analysis: First the energy 
densities of four little squares constituting a larger square are averaged. These averaged 
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energy densities are depicted in Fig. 9 as the second figure from the top of the left column. 
In other words, the original configuration has been smoothed or resolved on a rougher scale 
j=5. Some detail is obviously lost in this representation; no information about substructures 
of clumps living on the scale j = 5 can be deduced. This lost information can be recovered 
by keeping the difference between the resolutions j—5 and j=6; this difference information is 
illustrated in the picture to the right of the j—5 resolution picture. If this difference vanishes 
in some regions, no substructures are present; if on the other side those differences become 
large in other regions, then this indicates sizeable substructures on the scale j=6. 

These smoothing and differentiation operations are iterated through scales j=4, 3 down to 
2. In general, the difference between two successive smoothings of a configuration gives the 
information about the subclustering present at the involved scales. It is important to notice 
that the difference information on one scale is completely independent of (orthogonal to) 
the difference information on any other scale and therefore does not carry along redundant 
information. Also, no information is lost by keeping only the differences of smoothed repre- 
sentations. The sum of all difference representations together with the roughest smoothed 
configuration (j=2) recovers the original configuration belonging to the finest resolution 
scale. This is the two-dimensional analogue to Eq. (|3~9D . 

The smoothing operations at scale j are performed with the basis 

<t>jMk2( x ,y) = <\>{2~ 3 x - ki,2~ J y- k 2 ) = <fr(2~ 3 x - ki)(j){2- J y - k 2 ) = ^jkii^jkiiv) 5 (62) 



this is a direct product of two one-dimensional scaling functions given in eq. ( |3"o| ) as box 
functions. The difference between two successive smoothing operations can be completely 
expressed in terms of two-dimensional wavelets, namely 

^fLk 2 ^,y) = <f>jki(x)i/) jka (y) , 

V'Sxkafoy) = VW^^W?/) • (63) 



Here ipjk represents the one-dimensional Haar wavelet (p8|). Once again the dilation equations 
fl37|) govern the smoothing and differentiation transformations of the original configuration. - 
Referring to the previous paragraph, the difference information, i.e. the right column of Fig. 9, 
can be completely expanded in the two-dimensional Haar- wavelet basis (|63|). The correlations 
between the amplitudes of this multiresolution expansion are the wavelet correlations. Thus 
it becomes obvious that the wavelet correlations do characterize the subclustering of clumps. 

This interpretation works not only for the Haar wavelet basis. In Fig. 10 a two-dimensional 
wavelet analysis has been performed with the so-called Daubechies D12 wavelet [PH| . The 
decompositions fl62|) and fl63|) still hold; the only difference to the Haar-wavelet analysis is 
that 12 coefficients c m are used in the dilation equations ([37]). Again the left column of 
Fig. 10 shows density plots of a sequence of D12-smoothed approximations to the original 
configuration, whereas the middle column represents density plots of the wavelet transform, 
which form a sequence of mutually orthogonal details. In order to provide a better picture of 
the subclustering aspect of the wavelet transform, the details at various scales are exhibited 
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again in the right column, but with only two gray-values: black for regions where the detail 
function becomes negative, indicating local voids, and white for regions where the detail 
function is positive, signaling the appearance of clusters at the various scales. 

The borders between white and black regions in the right column of Fig. 10 are the zero- 
crossings of the wavelet transform. In signal analysis, these zero- crossings localize the signal 
sharp variation points at different scales and act as edge detectors |25j . The related wavelet 



transform maxima method EH] is used to compress and reconstruct signals efficiently 



4 Conclusions and outlook 

For the complex dynamics of multiparticle processes we study the fundamental differential 
correlation densities in order to extract as much important statistical information (scaling, 
clustering) as possible. This has been demonstrated with the help of some simple dis- 
crete hierarchical branching models. An evolution equation for the (multivariate) generating 
function has been formulated, from which the correlation densities follow recursively. An 
integration over the correlation densities yields the moments. 

Contrary to common practice, a careful conceptual distinction has to be made between 
forward (= theoretical) moments, which are obtained in evolution direction of the cascade 
process, and backward (= experimental) moments obtained in a backward analysis after the 
last cascade step has been performed. The former, often calculated from theory, can not be 
directly compared to the latter. Only the backward moments are accessible to experimental 
observation. We show that hierarchical branching processes with global density fluctuations 
on top of local density fluctuations do not show a rigorous scaling in the observed moments, 
although they are constructed by a purely selfsimilar iteration law. In other words, a devia- 
tion from scaling of the observed moments does not necessarily imply that the cascade process 
is not selfsimilar. This poses the provocative question: What is now the essence of a mo- 
ment (multifractal) analysis which only focuses on scaling indices? A deviation from perfect 
scaling provides important information on the cascade mechanism as we have demonstrated 
by various models. - In contrast to a moment analysis, the selfsimilar (forward) scaling can 
still be extracted from the differential correlation densities. 

However conventional correlation densities, especially higher orders, are hard to quantify. 
For their representation a compression is needed which reveals only the important (sta- 
tistical) information and removes redundancy. In this respect the wavelet transformation 
appears very appealing. Wavelets constitute a selfsimilar and orthogonal basis; in addition, 
due to their multiresolution properties, they dissect structures into details (clumps) living on 
different scales. As a consequence the wavelet transformed correlation densities of all orders 
become very sparse for hierarchically organized processes. The second order correlations 
diagonalize completely in the Haar-wavelet basis. Furthermore, for any order of the wavelet 
correlations, the diagonal elements show scaling exponents which can be related to the mul- 
tifractal dimensions. Beyond the multifractal analysis, those few off-diagonal elements of the 
higher order Haar-wavelet correlations which are nonzero and arise in bands show multiple 
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scaling and provide information on the subclustering structure of clumps such as correlations 
between small clusters living within larger clusters or voids. 

This paper provides a demonstration of the wavelet capabilities for correlation studies. 
However, the techniques presented here have to be further refined in order to become a pow- 
erful tool for experimental analysis. In this respect wavelet packets [O, |13[ might help to 



simplify the correlation structure even further. Also, in view of the recent success of corre- 
lation integrals P7[] , one should develop an analogue, so-called wavelet correlation integrals; 
this could set new standards of inferring information on clump correlations. Mathematically, 
there is the challenging problem to define the wavelet transformation on spheres which might 
have important applications in the context of angular intermittency. 

We believe that a true understanding of the intermittency phenomenon in multiparticle 
cascading is still a long way off. For general particle cascades it is not a priori clear that they 
follow a hierarchical evolution structure. Particle cascades should be formulated as point 
processes, which need not be hierarchically organized from the beginning. The relevant 
splitting functions should dictate if hierarchical structures could evolve dynamically. This 
might be the case for perturbative QCD. Work in this direction, including the use of the 
wavelet transformation, is presently being carried out. 

Apart from multiparticle physics, the above outlined wavelet technology might have ap- 
plications in other fields of physics. For example, wavelet (clump) correlations may shed 
further light on the long-standing question in intermittent fully developed turbulence of how 
eddies decay at small scales and would allow to gain new and supplementary insight for phe- 
nomenological modeling. In disordered solid state systems multifractality and localization 



of electron wavefunctions are discussed; a wavelet analysis appears to be suitable |28| . Last 
not least, the simplifying aspect of wavelet correlations might ignite the century-old inverse 
problem of how to extract dynamics from correlations with new fuel. 
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Figure Captions 



Figure 1: One possible tree configuration with random branching. At each cascade step 
all existing branches, which have evolved so far, might independently split into a left and a 
right subbranch with probability p or might not split with a probability 1 — p. 

Figure 2: Energy curdling for a binary cascade with random branching shown for three 
cascade steps. After the last cascade step (J = 3) each energy density is resolved at the 
finest scale. 

Figure 3: Reduced two-bin correlation density r klk2 = Pk 1 k 2 /Pk 1 Pk 2 °f the p-model (a), 
a-model (b), p-model with random branching (c) and the QCD-motivated cascade model 
(d). The splitting parameters and the number of cascade steps have been chosen to be 
a = (3 = 70 = 0.4 and J = 6, respectively; for case (c) the branching probability has been 
set to p — 0.5. 

Figure 4: Deviations from scaling of the backward moments M 2 (J,j) (upper two curves) 
and Haar wavelet correlations p^\\ ^ k \ (lower two curves) with respect to the scale index j 
for the p- (solid line) and a-model (dashed line). The number of cascade steps is J = 10 
and a = (3 = 0.4. 

Figure 5: Deviations from scaling of the normalized forward moments m 2 (j) (dashed line), 
the normalized backward moments M 2 (J,j) = M 2 (J,j)/(paP) 2 (solid line) and the normal- 
ized Haar wavelet correlations (jk) = p\jI) (jk)/ '{Po^Y (dotted line) with respect to the 
scale index j for the QCD-motivated cascade model. The number of cascade steps is J = 10 
and 70 = 0.4. 

Figure 6: Reduced Haar wavelet transformed two-bin correlation density vW k \ <j 2 k 2 ) = 

PUiki) {j 2 k 2 )i ^{Po^Y °f tli e p-model (a), a-model (b), p-model with random branching (c) and 
the QCD-motivated cascade model (d). The splitting parameters and the number of cascade 
steps have been chosen to be a = (3 = 70 = 0.4 and J = 6, respectively; for case (c) the 
branching probability has been set to p — 0.5. 

Figure 7: Deviations from scaling of M 2 (J,j) (upper curves) and P^k){jk) (l° wer curves) for 
the p-model with random branching with J = 10, a = 0.4, p — 1 (solid curve), p — 1/ (1 +a 2 ) 
(dashed curve) and p = 0.5 (dotted curve). 

Figure 8: The p- model evolution of the energy densities at three adjacent cascade steps. 
The marked arrows represent the difference in energy density of neighbouring bins and are 
equal to the corresponding Haar wavelet amplitudes 2ijk- The clump correlations relate the 
fluctuations of e.g. arrow 1 with those of arrow 2. 
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Figure 9: Multiresolution analysis with the Haar wavelet basis of one particular two- 
dimensional a-model realization at scale J = 6 down to j = 2. Left column: sequence 
of smoothing operations. Middle column: difference between two adjacent smoothed scales. 
Right column: clump structure of middle column emphasized by reduction to two gray values 
(white/black) for positive/negative regions in the difference information. 

Figure 10: Multiresolution analysis of 500 points from a Poisson transformed a-model 
realization in two dimensions with respect to the smooth compact Daubechies D12 wavelet. 
For further details see Fig. 9. 



31 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9503235vl 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9503235vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9503235vl 



This figure "fig2-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9503235vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9503235vl 



This figure "fig2-3.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9503235vl 



